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LONG-TERM  GOALS 

The  long-term  goal  of  this  research  is  to  develop  improved  parameterizations  of  internal  wave  driven 
mixing  in  global  circulation  models  where  such  small-scale  processes  are  not  resolved  (i.e.  occur  as 
subgrid-scale  processes). 

OBJECTIVES 

The  main  objectives  are:  (1)  to  perform  a  series  of  high-resolution  “numerical”  microstructure  profile 
studies  using  representative  field  scale  flows  in  the  coastal  ocean  to  investigate  the  relationship 
between  relevant  length  scales  and  time  scales  for  mixing  of  both  momentum  and  scalar  (density);  (2) 
to  formulate  parameterizations  for  diapycnal  (irreversible)  mixing  driven  by  internal  waves;  and  (3)  to 
explain  the  fate  of  nonliner  internal  waves  (NLIWs)  that  propagate  into  shallower  seas  as  a  result  of 
the  interaction  process  and  determine  whether  they  contribute  significantly  to  mixing  and  transport. 

APPROACH 

This  research  takes  a  take  multi-pronged  approach  that  involves  high-resolution  numerical  simulations, 
theoretical  analysis  and  model-data  comparisons.  The  first  major  thrust  of  this  work  is  to  perform 
highly  resolved  idealized  simulations  to  conduct  a  process-oriented  study  since  they  allow  for  specific 
internal  wave  induced  turbulent  mixing  processes  to  be  isolated  and  examined.  Furthermore,  they  will 
also  permit  the  exploration  of  a  broad  flow  and  environmental  parameter  space.  The  second  major 
thrust  of  this  research  project  is  to  perform  representative  field-scale  simulations  that  will  require 
subgrid-scale  modeling  to  resolve  the  turbulence.  Based  on  the  insights  gained  from  the  process- 
oriented  simulations  and  model-data  comparisons,  it  will  be  possible  to  evaluate  the  performance  of 
the  subgrid-scale  parameterizations  in  the  field-scale  simulations. 

WORK  COMPLETED 

For  FY2014,  we  have  focused  on  formulating  and  testing  a  simple  turbulence  parameterization  for 
stably  stratified  flows  in  the  proximity  of  boundaries  and  working  on  setting  up  and  running  process- 
oreinted  simulations.  We  have  also  focused  on  extending  the  findings  from  our  direct  numerical 
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simulations  (DNS)  of  homogeneous  stably  stratified  turbulence  (see  Mater  et  al.  2013)  to  field  scale 
observations.  Extension  of  this  published  work  is  now  in  progress  to  investigate  the  trends  in  ocean 
turbulence  data.  This  work  is  in  collaboration  with  Dr.  Lou  St.  Laurent  at  WHOI. 

RESULTS 

We  use  both  a  simple  zero-equation  turbulence  closure  scheme  as  well  as  the  standard  k-  s  model  to 
highlight  the  improvements  obtained  from  a  new  parameterization  for  the  turbulent  Prandtl  number 
that  we  have  derived.  The  turbulent  Prandtl  number  ( Prt )  is  the  linking  bridge  between  the  turbulent 
momentum  and  scalar  fluxes  in  the  context  of  Reynolds-averaged  Navier-Stokes  (RANS)  simulations. 
The  proposed  parameterization  for  Prt  takes  into  account  (for  the  first  time  to  our  knowledge)  the 
inhomogeneity  caused  by  the  wall  coupled  with  the  effects  of  density  stratification.  This  formulation  is 
given  in  equation  1  below  as: 


Pr,= 


1--L 

D 


+  Pr 

twd  0^  1  Lt  O’ 


(1) 


where  Rig  is  the  gradient  Richardson  number,  Rf  is  the  flux  Richardson  number,  z  is  height  measured 
from  the  bottom,  D  is  the  total  water  depth,  Prtwdo  is  the  difference  between  the  neutral  Prandtl  number 
at  the  wall  and  the  neutral  turbulent  Prandtl  number  for  a  homogeneous  shear  flow  ( Prt0 ). 

The  comparisons  (as  shown  in  figures  1  and  2)  of  the  one -dimensional  simulations  results  for  both  the 
mean  velocity  profile  and  the  mean  density  profile  with  direct  numerical  simulation  (DNS)  of  stably 
stratified  channel  flow  results  show  remarkable  agreement.  Also  shown  in  these  figures  are  the 
turbulent  Prandtl  number  formulations  proposed  by  Munk  and  Anderson  (1948)  and  a  homogeneous 
formulation  given  in  equation  (2)  below. 


?rt=-^  +  ?vt0.  (2) 

Rf 

Further  details  can  be  found  in  a  published  article  in  Journal  of  Geophysical  Research  (see  Karimpour 
and  Venayagamoorthy  2014a). 

The  next  phase  of  our  work  on  this  project  is  on  mixing  driven  by  breaking  internal  waves  using  highly 
resolved  numerical  simulations.  As  an  example,  Figure  3  shows  a  time  sequence  of  the  density  field 
(with  superimposed  velocity  vectors)  obtained  from  a  preliminary  high-resolution  laboratory  scale 
simulation  of  the  interaction  of  a  nonlinear  first-mode  internal  wave  interaction  with  a  steep  ridge 
topography.  The  rich  dynamics  of  the  interaction  process  is  evident  with  steepening  of  isopycnals  at 
the  leading  edge  of  the  wave  leading  to  local  kinematic  instabilities  that  cause  wave  breaking  as  well 
as  ejection  of  denser  fluid  over  the  ridge  crest.  At  the  rear  end,  the  formation  of  a  hydraulic  jump 
occurs  leading  to  enhanced  mixing  and  release  of  lee  waves.  These  initial  results  provide  some 
qualitative  insights  on  phenomena  associated  with  strong  wave  forcing  over  steep  topography  such  as 
the  double  ridge  system  in  Luzon  Strait  (e.g.  see  Alford  et  al.  2011).  Field  measurements  from  the 
ONR  DRI  on  Internal  Waves  in  Straits  Experiment  (IWISE)  cruises  in  2010  and  2011,  respectively 
indicate  overturning  length  scales  of  the  order  hundreds  of  meters  in  the  Luzon  Straits  at  certain 
locations  and  times  (e.g.  see  Alford  et  al.  2011  and  Louis  St.  Laurent  -  personal  communication).  It  is 
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not  clear  whether  such  large  overturns  are  indicative  of  actual  local  mixing  or  indicate  signatures  of 
dense  water  uplifted  higher  into  the  water  column  by  slope  convection  under  strong  wave  forcing.  Our 
ongoing  process-oriented  simulation  study  will  provide  insights  on  this  issue. 


Figure  1:  (left)  Velocity  profiles  and  (right)  density  profiles  for  Ri  r  =  60  obtained  from  the  zero- 
equation  closure  scheme  using  (a,  b)  the  Prt  formulations  given  by  equation  (2)  and  the  formulation 
of  Munk  and  Anderson  (1948)  and  (c,  d)  the  modified  Prt  formulation  (equation  ( 1 ))  compared  with 
the  channel  flow  DNS  data  of  Garcia-Villalba  and  del_Alamo  (2011). 

IMPACT/APPLICATIONS 

The  relevance  and  broader  scientific  impacts  of  this  research  project  come  from  the  importance  of 
turbulent  mixing  and  transport  processes  in  geophysical  and  environmental  flows.  The  expected  results 
from  this  research  will  enhance  our  basic  understanding  of  turbulent  mixing  and  transport  in  the 
coastal  ocean  induced  by  internal  waves  interacting  with  topography.  This  topic  is  a  very  active  area  of 
research  in  the  oceanic  research  community.  A  major  expected  outcome  of  this  research  is  the 
formulation  of  simple  but  yet  accurate  parameterizations  of  mixing  and  transport  that  will  be  valuable 
in  developing  better  subgrid  scale  models  for  use  in  large  scale  numerical  circulation  models  of  the 
coastal  ocean. 
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Figure  2:  (left)  Velocity  profiles  and  (right)  density  profiles  for  Rir  =  60  obtained  from  the  standard 
k-e  closure  scheme  using  (a,  b)  the  Prt  formulations  given  by  equation  (2)  and  the  formulation  of 
Munk  and  Anderson  (1948)  and  (c,  d)  the  modified  Prt  formulation  (equation  ( 1 ))  compared  with 
the  channel  flow  DNS  data  of  Garcia-Villalba  and  del_Alamo  (2011). 

RELATED  PROJECTS 

The  PI  has  another  ONR  funded  project  on  improved  turbulence  parameterizations  for  oceanic  flows 
where  the  goal  is  to  develop,  implement  and  test  different  turbulence  parameterizations  in  numerical 
models.  The  emphasis  of  this  work  is  on  developing  robust  parameterizations  that  are  applicable  for  a 
wide  range  of  oceanic  flow  conditions.  Hence  there  is  some  natural  overlap  between  these  two 
projects.  Also,  the  PI  (through  a  related  topic  funded  partly  through  an  NSF  grant)  has  been  able  to 
investigate  boundary-layer  turbulence  especially  very  close  to  the  wall.  Our  findings  have  been 
recently  published  in  two  articles  in  the  Journal  of  Fluid  Mechanics  (see  Karimpour  and 
Venayagamoorthy  2014b,  2013).  We  have  extended  this  work  to  wall-bounded  stratified  flows  and 
therefore  the  efforts  are  synergistic  to  this  project. 
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Figure  3:  Time  sequence  of  density  field  and  velocity  vectors  from  preliminary  idealized  numerical 
simulations  to  illustrate  the  rich  flow  dynamics  that  occurs  when  a  low-mode  internal  wave 
impinges  on  a  steep  ridge.  Note  how  dense  fluid  is  advected  up  the  ridge  during  the  flood  phase  of 
the  tide  (wave  forced  from  left  end  boundary  which  is  truncated  for  clarity)  and  the  generation  of  an 
internal  hydraulic  jump  and  the  lee  wave  release  during  the  ebb  (slack)  phase  of  the  wave. 
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